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Abstract 

Stochastic differential equations are often simulated with the Monte Carlo Euler method. Conver- 
gence of this method is well understood in the case of globally Lipschitz continuous coefficients of the 
stochastic differential equation. The important case of superlinearly growing coefficients, however, 
has remained an open question. The main difficulty is that numerically weak convergence fails to 
hold in many cases of superlinearly growing coefficients. In this paper we overcome this difficulty and 
establish convergence of the Monte Carlo Euler method for a large class of one-dimensional stochastic 



differential equations whose drift functions have at most polynomial growth. 

B 

1 Introduction 

(N 

^ Many applications require the numerical approximation of moments or expectations of other functionals 

of the solution of a stochastic differential equation (SDE) vifhose coefficients are superlinearly growing. 



Moments arc often approximated by discretizing time using the stochastic Euler scheme (see e.g. 11 



15 , 21 ) (a.k.a. Euler-Maruyama scheme) and by approximating expectations with the Monte Carlo 
method. This Monte Carlo Euler method has been shown to converge in the case of globally Lipschitz 
continuous coefficients of the SDE (see e.g. Section 14.1 in jllj and Section 12 in [15| ). The important 
case of superlinearly growing coefficients, however, has remained an open problem. The main difficulty 
is that numerically weak convergence fails to hold in many cases of superlinearly growing coefficients; 
j>! see [t]. In this paper we overcome this difficulty and establish convergence of the Monte Carlo Euler 

• ^ method for a large class of one-dimensional SDEs with at most polynomial growing drift functions and 

rN with globally Lipschitz continuous diffusion functions; see Section [2] for the exact statement. 

For clarity of exposition, we concentrate in this introductory section on the following prominent 
example. Let T e (0, oo) be fixed and let {Xt)te[Q.T] be the unique strong solution of the one-dimensional 
SDE 

dXt = -Xf dt + a dWt, Xo = xq (1) 

for all t G [0,T], where (W^t)te[o.T] is a one-dimensional standard Brownian motion with continuous 
sample paths and where a € [0,00) and a;o G M are given constants. Our goal is then to solve the 
cubature approximation problem of the SDE ([T]). More formally, we want to compute moments and, 
more generally, the deterministic real number 

e[/(^t)] (2) 
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for a given smooth function / : E — M whose derivatives have at most polynomial growth. 

A frequently used scheme for solving this problem is the Monte Carlo Euler method. In this method, 
time is discretized through the stochastic Euler scheme and expectations are approximated by the Monte 



Carlo method. More formally, the Euler approximation (Y^), 
SDE ([T]) is defined recursively through = xq and 



e{o,i....,Af} of the solution {Xt)f^rQ j^i of the 



Y 



N 



T 

N 



Y, 



(3) 



for every n e {0, 1, . . . , iV - 1} and every e N := {1, 2, . . .}. Moreover, let F„^'", n e {0, 1, . . . , N}, 
e N, for m S N be independent copies of the Euler approximation defined in (|3|. The Monte Carlo 
Euler approximation of ^ with N & N time steps and Af € N Monte Carlo runs is then the random real 
number 

1 / *^ \ 



M 



E 



f{y, 



N 



(4) 



In order to balance the error due to the Euler method and the error due to the Monte Carlo method, it 
turns out to be optimal to have M increasing at the order of N"^; see [2]. We say that the Monte Carlo 
Euler method converges if 



lim 

JV-yoc 



E 



1 



E 

m— 1 



N,'i 
N 







(5) 



whose derivatives have at most polynomial 



holds almost surely for every smooth function / : M - 
growth (see also Appendix A.l in [4]). 

In the literature, convergence of the Monte Carlo Euler method is usually established by estimating 
the bias and by estimating the statistical error (see e.g. Section 3.2 in [20| ). More formally, the triangle 
inequality yields 



E 



1 



N.ni\ 

N ' ) 



< 



E 



-E 



f{Yi 



N J 



approximation error of the 
Monte Carlo Euler method 



absolute value 
of the bias 



E 



.f{Yj 



N J 



1 



7V2 



f{Yi 



N,m\ 
N ) 



(6) 



statistical error 



for every A^ e N and every smooth function / : M — > M whose derivatives have at most polynomial growth. 
The first summand on the right-hand side of Q is the absolute value of the bias due to approximating the 
exact solution with Euler's method. The second summand on the right-hand side of ([6]) is the statistical 
error which is due to approximating an expectation with the arithmetic average over independent copies. 
The bias is usually the more difhcult part to estimate. This is why the concept of numerically weak 
convergence, which concentrates on that error part, has been studied intensively in the literature (see 
for instance [li], [l2], [l5], (20], [l], [s], [o], [l^ or Part VI in [IT]). To give a definition, we say that 
the stochastic Euler scheme converges in the numerically weak sense (not to be confused with stochastic 
weak convergence) if the bias of the Monte Carlo Euler method converges to zero, i.e., if 



lim 



E 



E 



![Yi 







(7) 



holds for every smooth function /: M — > M whose derivatives have at most polynomial growth. If the 
coefficients of the SDE are globally Lipschitz continuous, then numerically weak convergence of Euler's 
method and convergence of the Monte Carlo Euler method is well-established; see e.g. Theorem 14.1.5 
and Section 12 in 



11 
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The case of superlinearly growing coefficients is more subtle. The main difficulty in that case is that 
numerically weak convergence usually fails to hold; see [7j for a large class of examples. In particular, 
the sequence E[(Y^)^] , A^ S N, of second moments of the Euler approximations ^ diverges to infinity 
if CT > although the second moment E[(ArT)^] of the exact solution of the SDE (fl| is finite and, hence, 
we have 



lim 

iV->c» 



E 



{XtY 



E 



(8) 



2 



instead of ([t]) . The absolute value of the bias thus diverges to infinity in case of SDEs with superlinearly 
growing coefficients. This in turn implies divergence of the Monte Carlo Euler method in the mean square 
sense, i.e., 



E 



E 



1 

iV2 



E 



[Xt? 



E 



N 



Var 



squared bias — >-oo 




— )■ OO 



mean square error of the 
Monte Carlo Euler method 



varianec of the Monte 
Carlo Euler method 



(9) 

as iV — ^ oo. Clearly, the mean square divergence ([9]) does not exclude the almost sure convergence ([5|. 
Indeed, the main result of this article proves the convergence ^ of the Monte Carlo Euler method. For 
proving this result, we first need to understand why Euler's method does not converge in the sense of 
numerically weak convergence. In the deterministic case, that is, ([I]) and ([s]) with ct = 0, the Euler 
approximation diverges if the starting point is sufficiently large. This divergence has been estimated 
in [t] and turns out to be at least double-exponentially fast. Now in the presence of noise ((j > 0), the 
Brownian motion has an exponentially small chance to push the Euler approximation outside of [— A^, N]. 
On this event, the Euler approximation grows at least double-exponentially fast due to the deterministic 
dynamics. Consequently, as being double-exponentially large over-compensates that the event has an 
exponentially small probability, the L^-norm of the Euler approximation diverges to infinity and, hence, 
numerically weak convergence fails to hold. 

Now we indicate for example ([ij with — Q and d — 1 why the Monte Carlo Euler method converges 
although the stochastic Euler scheme fails to converge in the sense of numerically weak convergence. 
Consider the event fijv := {supo<t<7i |Wt| < yjN/{2T)} and note that the probability of {^nY is 
exponentially small in iV g N. The key step in our proof is to show that the Euler approximation does 
not diverge on Vl^ as iV € N goes to infinity. More precisely, one can show that the Euler approximations 
([3]) are uniformly dominated on Vljq by twice the supremum of the Brownian motion, i.e., 



sup ( la„ Ir^l 



< 2 ( sup \Wt 

yO<t<T 



holds. Consequently, the restricted absolute moments are uniformly bounded 



sup E 

JVGN 



N\P 



< 2P -E 

for all p e [1,00). This estimate complements the divergence 



sup \Wt 

a<t<T 



< 00 



lim E 

N—^oo 



Y, 



N\P 
N 



00 



(10) 



(11) 



(12) 



for all p £ [1,00), which has been established in [t]. Now once the restricted absolute moments are 
uniformly bounded, an adaptation of the arguments of the globally Lipschitz case leads to the modified 
numerically weak convergence 



lim E 

AT— >oo 



,f{Yi 



N ) 



= E 



(13) 



for every smooth function / : M — > M whose derivatives have at most polynomial growth, see Lemma |4.6| 
By substituting this into an inequality analogous to ^ and by using the exponential decay of the 
probability of {flNY in N G N, one can establish convergence of the Monte Carlo Euler method. Note 
that a domination as strong as ( 10 1 holds for more general non-increasing drift functions if the diffusion 
function is identically equal to 1. For more general drift and diffusion functions, however, both flj^r and the 
dominating process are more complicated in that they depend on the Euler approximation. Nevertheless, 
the dominating process can be shown to have uniformly bounded absolute moments; see Section |4] for 
the details. 

Our main result, Theorem |2.1| below, establishes convergence of the Monte Carlo Euler method for 
SDEs with globally one-sided Lipschitz continuous drift functions and with globally Lipschitz continuous 
diffusion functions. Moreover, the coefficients of the SDE are assumed to have continuous fourth deriva- 
tives with at most polynomial growth, see Section |2] for the exact statement. The order of convergence 
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turns out to be as in the globally Lipschitz case. In that case, the stochastic Euler scheme converges 
in the sense of numerically weak convergence with order 1. The Monte Carlo simulation of E[/(Y^)] 
with M independent Euler approximations has convergence order ^ — . For a real number r > 0, we write 
r— for the convergence order if the convergence order is better than r — e for every arbitrarily small 
e £ (0, r). We therefore choose M — N"^ in order to balance the error arising from Euler's approximation 
and the error arising from the Monte Carlo approximation. Both error terms are then bounded by a 
random multiple of A^('^~i) with e € (0, 1). Since 0{M ■ N) = 0{N^^ function evaluations, arithmetical 
operations and random variables are needed to compute the Monte Carlo Euler approximation Q, the 
Monte Carlo Euler method converges with order 



of global Lipschitz coefhcients of the SDE (see 21). Theorem 



with res pect to the computational effort in the case 
shows that |— is also the convergence 



2.1 



SDE. Simulations support this result, see 



order in the case of superlinearly growing coefficients of the 
Section [3l 

Let us reconsider the standard splitting ^ of the approximation error into bias and statistical error. 
Theorem 2.1 of |7 implies that the absolute value of the bias diverges to infinity as N ^ oo. This together 
with our Theorem 2.1 below yields that also the statistical error diverges to infinity. More formally, we 

see that 



E 



1 
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m— 1 
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m — 1 



N.; 
N 



(14) 



statistical error — >oo 



P-a.s. as — >■ oo for every smooth function / : M — > M with at most polynomially growing derivatives 
and with f{x) > c\x\'^ — 1/c for all x e M and some c G (0, oo). This emphasizes that the standard 
splitting of the approximation error of the Monte Carlo Euler method into bias and statistical error is 
not appropriate in case of SDEs with superlinearly growing coefficients. 



2 Main result 

We establish convergence of the Monte Carlo Euler method for more general one-dimensional diffusions 
than our introductory example ([I]). More precisely, we pose the following assumptions on the coefficients. 
The drift function is assumed to be globally one-sided Lipschitz continuous and the diffusion function 
is assumed to be globally Lipschitz continuous. Additionally, both the drift function and the diffusion 
function are assumed to have a continuous fourth derivative with at most polynomial growth. 

We introduce further notation for the formulation of our main result. Fix T G (0, oo) and let (il, J^, P) 
be a probability space with a normal filtration {J^t)t<^[o t] ■ ■ — >■ K, m G N, be a sequence of inde- 

pendent, identically distributed Jx)/'S(M)-measurable mappings with E[|^'^^^|p] < oo for every p € [1, oo) 
and let : [0, T] x — > E, to G N, be a sequence of independent scalar standard iJ^t)t,= [o T^j-Brownian 

motions with continuous sample paths. Furthermore, let /i, a : K — >■ M be four times continuously differ- 
entiable functions. Generalizing ([T]), let {Xt)te[o.T] be a one-dimensional diffusion with drift function fi 
and diffusion function a. More precisely, let X : [0, T] x $7 — >■ K be an (up to indistinguishability unique) 
adapted stochastic process with continuous sample paths which satisfies 

Xt^C^''^+ I ^i{Xs)ds+ f a{X,)dW^'^ 



1 (15) 

for every t £ [0,r]. The functions fi and are the infinitesimal mean and the infinitesimal variance 
respectively. 

Next we introduce independent versions of the Euler approximation. Define J^/,B(M)-measurable 
mappings Y^'"" : f7 ^ M, n e {0, 1, . . . , TV}, TV e N, to e N, by Yo^'"\lu) C^'")(a;) and by 

y„^rM Y^^^^iu^) + ^ ■ M(i;^''"(a.)) + a(r„^'™(a.)) • (w<l^{u^) - W^^\co)^ (16) 

for every n G {0, 1, . . . , N— 1}, N ^ N and every m S N. Now we formulate the main result of this article. 
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Theorem 2.1. Suppose that f, ^, a: 



are four times continuously differentiable functions with 



(x) + cr(")(j;) <L(l + |a;|'') Vx e 



(17) 



for every n G {0, 1, • • . , 4}, where L G (0, oo) and 6 G (1, oo) are fixed constants. Moreover, assume that 
the drift coefficient is globally one-sided Lipschitz continuous 



{x - y) ■ {p{x) - ^i{y)) < L {x - y) Va;,ye 
and that the diffusion coefficient is globally Lipschitz continuous 

Hx) - (j{y)\ < L\x ~ y\ Va;,yeM. 



(18) 



(19) 



Then there are T /B{[Q,oo)) -measurable mappings : — > [0, oo), e £ (0,1), and a set ^ £ J- with 
F[Cl] = 1 such that 



E 



L ^ /(f#-(^)) 



iV2 



1 



Ar(i-£) 



(20) 



holds for every w G fJ, G N and every e G (0, 1). 



The proof is deferred to Section |4j For further numerical approximation results for SDEs with super- 
linearly growing coefficients, see e.g. [Gj, [13] , 16 , 19 and the references in the introductory section 
of[7). 

Note that the assumptions of Theorem |2.1| ensure the existence of an adapted stochastic process 



X : [0, r] X — M with continuous sample paths which satisfies ( 15 ) and 



E 



sup \Xt\^ 

0<t<T 



< OO 



(21) 



2.1 



for all p G [l,oo) (see Theorem 2.6.4 in fsj). Therefore, the expression E[/(Xt)] in (201 in Theorem 
is well-defined. 

Since 0{N^^ function evaluations, arithmetical operations and random variables are needed to com- 
pute the expression (j2m=i /(^jv''™('^))) f^O) for lo £ fl, Theorem |2 . 1 1 shows that the Monte Carlo 

Euler method converges under the above assumptions with order with respect to the computational 
effort. This is the standard convergence order as in the global Lipschitz case (see e.g. [2]). 



3 Simulations 

In this section, we simulate the second moment of two stochastic differential equations. First we simulate 
the stochastic Ginzburg-Landau equation with multiplicative noise, which we choose as there exists an 
explicit solution for this SDE. Let {Xt)te[o,i] be the solution of 



dXt ^{\xt- X't ]dt + Xt dWu Xo = 1 



(22) 



for all t G [0, 1]. The exact solution at time 1 is known explicitly (see e.g. Section 4.4 in [11| ) and is given 
by 

exp(W^i) (23) 



Xi = 



1 + 2 J^cxp{2Ws)ds 



The exact value of the second moment E[(Xt)^] ~ E[(Xi)^] is not known. Instead we use the exact 
solution ( 23 ) at time 1 to approximate the second moment. For this, we approximate the Lebesgue integral 
in the denominator of ( 23 ) with a Riemann sum with 3 • 10'^ summands. Moreover, we approximate the 



second moment at time 1 by a Monte Carlo simulation with 10^ independent approximations of Xi. This 
results in the approximate value E[(Xi)^] w 0.4945. 
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TV = 2" 


TV = 21 


iV = 2^ 


= 2^ 


A^ = 2* 


1.1379 


0.9118 


0.4258 


0.2942 


0.4386 


iV = 2^ 


= 2^ 


N = 2'^ 


= 2^ 


A^ = 2" 


0.4641 


0.4663 


0.4859 


0.4904 


0.4935 



Table 1: Monte Carlo Eulcr approximations |24]l of E [{Xif] of the SDE ([22| 



Next we approximate the second moment at time 1 with the Monte Carlo Euler method. We will 
sample one random a; e and calculate the Monte Carlo Euler approximations for this w e for different 
discretization step sizes. More precisely, Table [l] shows the Monte Carlo Euler approximation 



1 ^ 

— y f 



-trN. 771/ \ 



(24) 



of the second moment at time 1 of the SDE (22 1 for every N £ {2^, 2^, 2^, . . . , 2^} and one random cj G f^. 
In Figure [l| the approximation error of these Monte Carlo Euler approximations, i.e., the quantity 



ApproKimation error of E(X^]^ of the stochastic Ginz burg -Landau equatiorifor N^^*^,?', .,2^ 

'~ ' Appro>;imati on error 

"O" Order lines 1/6, 1/3 and 1/2 




10° 10' 1Cf 10^ 10" lO'" 10'' 

Compulational efforl 



Fi gure 1 '. Approximation error pSl l of the Monte Carlo Eulcr approximations ([24]| of E[(Xi)2] of the SDE ([22} 



0.4945 



-xrN.rn ( \ 



(25) 



is plotted against A"^ for every A" € {2°, 2^, 2^, . . . , 2^}. Note that is the computational effort up 
to a constant. The three order lines in Figure [l] correspond to the convergence orders | and \. 
Hence, Figure [T] indicates that the Monte Carlo Euler method converges in the case of the stochastic 
Ginzburg-Landau equation (22 1 with its theoretically predicted order | — . 

Next we simulate our introductory example. Let (Xi)(£[o.T] be the solution of the SDE ([T]) with 
T = 1, CT = 1 and = 0. The SDE Q thus reads as 



dXt = -Xl dt + dWt, Xo^O 



(26) 



for all t e [0,1]. Here there exists no explicit expression for the solution or its second moments. As 
an approximation of the exact value E[(Ari)^] , wc now take a Monte Carlo Euler approximation with a 
larger A^. We choose A^ = 2^^ and obtain the value 0.4529 ~ E[(Xi)^] as an approximation of E[(Xi)^] . 
Table [2] shows the value of the Monte Carlo Euler approximation 
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TV = 21 


iV = 2^ 


= 2^ 


A^ = 2* 


1.4516 


0.5166 


0.4329 


0.5308 


0.4285 


N = 2^ 


= 2^ 


iV = 27 


A^ = 2^ 


A^ = 2" 


0.4452 


0.4602 


0.4517 


0.4548 


0.4537 



Table 2: Monte Carlo Euler approximations ^ of E[(Xi)2] of the SDE (|26| 



(27) 



of the second moment at time 1 of the SDE (26 1 for every N ^ {2°, 2^, 2^, . . . , 2^} and one random w G fJ. 
In Figure [2j the approximation error of these Monte Carlo Euler approximations, i.e., the quantity 



Approximation error of E(X )^ of equatiori [1J for N =2°, 2' 2^ 




1[f 10' 10' 10' 

Compiiatiorial effort 



10 10 10 



Figure 2: Approximation error ([28]l of the Monte Carlo Euler approximations ([27|l of E[(Xi)2] of the SDE ([26j 



1 



(28) 



is plotted against A"^ for every A' G {2", 2^, 2^, . . . , 2^}. Note that N'^ is the computational effort up to a 
constant. The three order lines in Figure [2] correspond to the convergence orders g, \ and \. T herefore, 
Figure [2] suggests that the Monte Carlo Euler method converges in the case of the SDE ( 26 1 with its 
theoretically predicted order | — . 



4 Proof of Theorem 2.1 



First we introduce more notation. Recall the standard Brownian motion W^^^ : [0, T] x R and the 
Euler approximations Yn'^ : 17 — M, n G {0, 1, . . . , A^}, A^ G N, from Section[2j Throughout this section, 
we use the stochastic process : [0, T] x O E and the J^/S(M)-measurable mappings Y^ : -> K, 
n G {0, 1, . . . , A^}, A^ G N, given by 



Wt{^) := W[^\^) and Y^^ {u:) r„^-i(..) 
for every t G [0, T], n G {0, 1, . . . , A^}, lu G and every A^ G N. 



(29) 
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4.1 Outline 



For general drift and diffusion functions, our proof of Theorem |2 . 1 1 somewhat buries the main new ideas. 
To explain these ideas, we give now a very rough outline (the precise estimates and assertions can be 
found in Subsections 4.2|4.7 below). The main step will be to establish uniform boundedness of the 
restricted absolute moments of the Euler approximations 



sup E 



\VN\ 



< oo 



(30) 



where (r2jv)7veN is a sequence of events whose probabilities converge to 1 sufficiently fast. From here, 
one can then adapt the arguments of the global Lipschitz case to derive the modified numerically weak 



convergence (13) and to obtain Theorem 2.1 



The idea behind (30) is now explained on the example of negative cubic drift and multiplicative noise. 
Formally, we consider ^{x) — -x^, <t{x) = X for aU X e M and ^^^^{oj) = 1 for all oj e n. The Euler 
approximation ( 29 ) is then given by Y^^ = 1 and 



N 



N 



T 
N 



N 



(31) 



for all n e {0, 1, 

(Xk)k£{o.i,...,N} does not change sign until and including the n-th approximation step for some n G 



— 1} and all N G N. Fix e N and assume that the Euler approximation 
n 

> for all fee {0, 1, ... , n}. Then, using 1 + x < exp(a;) 



{0, 1, . . . , N} which we fix for now, that is, Yj^ 
for all a; S M, we have 



_ _ ^ (vN \3 



Y, 



N 

k-1 



.( 

1 (l + WkT - Wik-Dr) 

V " N / 

< Y^_j^exp(WkT - W 



(32) 



and iterating this inequality shows 



Yk < yk-2 expfT4^ (fc-i)T - W ik-2)A expf VK« 
Y^_^ expfw^fcT - W ik-2)A 
<...<y(f exp(M^^) 



W {k-l)T 



(33) 



for all k e {0, 1, . . . , n]. Thus the Euler approximation (i^jf )feg{o,i,....n} is bounded above by the domi- 
nating process (-D^)j.g{o,i,...,n}- This dominating process has uniformly bounded absolute moments. So 
the absolute moments of the Euler approximation can only be unbounded if the Euler approximation 
changes its sign. Now if Yk happens to be very large for one fc S {0, 1, . . . , — 1}, then Yk+i is negative 
with absolute value being very large because of the negative cubic drift. Through a sequence of changes 
in sign, it could happen that the absolute value of the Euler approximation increases more and more. To 
avoid this, we restrict the Euler approximation to an event 17 jv on which the drift alone cannot change the 
sign of the Euler approximation. On 0,^, the Euler approximation changes sign only due to the diffusive 
part. As the diffusion function is at most linearly growing, these changes of sign can be controlled. In 
between consecutive changes of sign, the Euler approximation is again bounded by a dominating process 
as above. Through this pathwise comparison with a dominating process, we will establish the uniform 
boundedness (30) of the restricted absolute moments. For the details, we refer to Lemma 4.2 which is 



the key result in our proof of Theorem |2.1 



4.2 Notation and auxiliary lemmas 



2.1 



In order to show Theorem 

n e {0, 1, . . . , iV} and everylVe N and let the B{M.) IB{. 



the following objects are needed. First of all, define 

measurable mapping it : M — > 



a{x) 



{a{x)~a(0)) 



x^O 
a; = 



: ^ for every 
be given by 



(34) 
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for all x e M. Moreover, let the J^/S(M)-measurable mappings a^'' 
by 

TL 



and 



M and P^^-'' 



be given 
(35) 

(36) 



1} and every N,m e N. For simplicity we also use , : 51 



for every uj £ fl, n € {0, 1, . . . , iV 
given by {oj) := a^'^{ui) and (w) := P^'^{uj) for every w G fi, n € {0, 1, . . . , iV — 1} and every 
N G N. Using these ingredients, we now define the dominating process. Let the J-"/S(M)-measurable 
mapping : — > M be given by 



D 



N.m 



(^) := e(^^- (r|/i(0)| + |a(0)| 

,n(yf-(c.))e(5:n.V."f'' 



w — 1 
k—v 



e('"'(w)| + 

("))/3f''"(c^) 



(37) 



for every lu £ fl, v,w E {0,1, ... , N} and every N,m £ N, where sgn: M — > {—1, 1} is given by sgn(a;) := 1 
for every x £ [0,oo) and sgn(x) := —1 for every x £ (— oo,0). As usual, X^ztLt;^ ~ for every 
V £ {w,w + 1,...,N}, w £ {0,1,..., N}, iu £ n and every N,m £ N. Note that D^]^ : n R only 
depends on the Brownian motion : [0, T] x 57 — >• M and the initial random variable 5^™-' : fl — ^ 

M for every v,w £ {0,1,..., A^} and every N,m £ N. Therefore, D^'^, m g N, is a sequence of 
independent random variables for every v,w £ {0, 1, ... , N{ and every A € N. We will show that the 
Euler approximation is dominated by the dominating process since the last change of sign. More formally, 
let the F /V{{0, 1, . . . , n})-measurable mapping : 51 {0, 1, . . . , n} be given by 



[uj) max [ {0} \J Ik £ {1,2, ... ,n} 



sgn(r,^Ii(^))^sgn(yf (..)) 



(38) 



for every uj £ VL, n £ {0,1, ... , N{ and every A S N. The random time r„ : 51 — )■ {0, 1, . . . , n} is the last 



4.2 



time of a change of sign of , k £ {0, 1, . . . , n}, for every n £ {0, 1, . . . , A^} and every N £f^. Lemma 
below shows that \Y^\ is bounded by D^n^^ on a certain event 5lAr^„ for every n £ {0, 1, ... , A^} and 
every A^ e N. Next we define these events i^N.n, n £ {0, 1, . . . ,N{, A^ G N, such that the drift alone 
cannot cause a change of sign and such that the increment of the Brownian motion is not extraordinary 
large. Let the real number rjv £ [0, oo) be given by 



rN 



Ai 

IT' 



max 0, 



A' 



r(sup,g[-i4]lM'(-s)l+3i 



for every N £'H. Now define sets ^N,m ^^i^N G -7^ by 



= < w e 51 




sup 

v,w£{0,l,. 



sup 
fce{o,i,---,i 



-1} 



WfN (w) 



(w) 



sup \D^'^{uj)\ < tn, sup 

v,w£{Q,l,...,N} ne{0,l,...,iV-l} 



< N- 



57 AT :— 51jv^Ar — 

for every n £ {0, 1, ... , A^} and every N,m £ N. Finally, define 51 G by 

/ oo 

^ - u n n 



(39) 

(40) 

(41) 
(42) 



,NGN M=N m=l 



(al 


uj£n 







sup 

NGN 



Ef^i(in^M-/(y^'"(^))-E[in./(y^)] 



(43) 



< oo 
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Note that (l is indeed in J^. Moreover, we write H^H^p (E[|Z|^])'' E [0, oo] for all p S [l,oo) and all 
J^/S(M)-measurable mappings Z: — M. Our proof of Theorem 2.1 uses the following lemmas. 

Lemma 4.1 (Burkholder-Davis-Gundy inequality). Let N E N and let Zi, . . . , Zm ■ ^ 

measurable mappings with E|Z„p < oo for all n E {1, . . . , N} and with E [Z„+i|Zi, . . 
n € {1, . . . , iV — 1}. Then we obtain 



■ R be J^/B{R)- 
Zn] = for all 



||Zi + . . . + Z^h, <Kp- (\\Zi\\% + ... + \\Zn\\1,) ■ 

for every p G [2,oo), where Kp, p G [2,oo), are universal constants. 
The following lemma is the key result in our proof of Theorem |2.1| 



Lemma 4.2 (Dominator Lemma). Let F„ : 



D 



N.l 



n 



(44) 



and fijv,n G J' for 



n,me {0, 1, . . . ,iV} and N e N be given by (29), (37 1, (381 and (40). Then we have 



\Y'^{lo) \ < D% , (w) 
, iV} and every iV £ N. 



for every lo £ ^N.n, n G {0, 1, 



(45) 



The domination ( 45 ) might not look helpful at first view since ^ depends on the Euler approximation 

and since is in general not a stopping time for n e {1, . . . , A^} and e N. However, the dependence of 
the dominating process on the Euler approximation can be controlled as a is bounded and the dependence 



of on for all n S {0, 1, . . . , iV} and all iV e N is no problem as -D„ J can be controlled uniformly 

in u, w e {0, 1; . . . , A^} and G N. This is subject of the following lemma. 

Lemma 4.3 (Uniformly bounded absolute moments of the dominator). Let D^^: 



n,m (z {0, 1, . . . , A'^} and N E N be given by (37). Then we have 



sup E 

N<£N 



N,1\P 



v,we{0.1,....N} 



< OO 



[0,oo) for 



(46) 



for all p G [1, oo). 



From Lemma [4.2| and from Lemma [4.3[ we immediately conclude that the restricted absolute moments 
of the Euler approximation are uniformly bounded. 

Corollary 4.4 (Bounded moments of the Euler approximation). Let the Euler approximation : — > 



for n G {0, 1, . . . , A^} and N €z N be given by (29). Then we have 



sup sup E 

Nmne{0,l,...,N} 

for every p £ [1, oo). 

Next we estimate the probability of rijv for every large A^ G N. 



< oo 



(47) 



Lemma 4.5 (Full probability). Let 51^? G J' for N E N and Q. E F be given by (42| and (43). Then we 
have that 



sup 



N) 



< oo 



and 



= 1 



(48) 



holds. 



Using the above lemmas, we establish the following modification of numerical weak convergence. 



Lemma 4.6 (Modified weak convergence). Let X : [0,T] x O 



Oat G J" and : n 



for 



G {0, 1, . . . , AT} and N e N be given by (15), (42| and (29). Then we obtain that 



sup ( N 

NeN 



E 



E 



< oo 



(49) 



holds. 



While the proof of Lemma 4.1 can be found in Theorem 6.3.10 in fl8l, the proofs of Lemma 4.2 



Lemma 4.3 Lemma 4.5 and Lemma 4.6 are given in Sections 4.3|4.7 below 
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4.3 Proof of Theorem [2J] 



Proof of Theorem 2.1 Consider the J-/B ([0, oo))-measurable mapping Z: il — >■ [0, oo) given by 



:= sup 
Nefi I 



(UJ) 



J2 \fiY^^"\^)) 



1 - 1 



(^) 



(50) 



for every uj £ Q. Finiteness in (50 1 follows from 



lim (l/„oo „m2 o^^U'^) ) = !(', 1 o^U^) ^ = 1 (51) 

for every w G Jl. Moreover, define J^/S([0, oo))-measurable mappings i?e : — > [0,oo), e e (0, 1), by 



i?e(w) := sup 



Ar(i+e) 



(52) 



for every a; G and every e £ (0, 1). By definition of ft, the mappings i?^ : — >■ [0, oo), e G (0, 1), are 
also finite. Additionally, let the real number C £ [0, oo) be given by 



C := sup ( TV 

NeN 



E 



E 



lf2N • fiY^] 



sup ( • P 



(53) 



which is finite due to Lemma |4.5| and Lemma |4.6[ Moreover, we have 



E 



< 



iV2 



E 



E 



E 



E 



E 



iV2 



for every cj e and every N £N. Using (|53|), we then obtain 

fiXr 



E 



< 



/ w 



E 



1 



L m — 1 

at" 



1 

N 



rn—l 



E if{Y^'"'i^))-ln- (-)-E 
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and 



E 



< 



1 / ^ 



E 



Ar(i+e) 



m— 1 



E 



sup (l-lo™(a;)) 

me{l,2,...,Af2} 



for every € fJ, A'^ G N and every e G (0, 1). Hence, we get 



E 



fiXr) 



1 / ^ 



< E 



TV 



^(ll/(-^-"<^'))l)(. 



sup (l-ln-(w)) 
ie{i,2,...,A'2} 



and 



E 



< E 



m—1 



E 



1/(^2 



1 / ^ , 

ME /(^#'™( 



Ar2 



1- inf lo-M 



(54) 



(55) 



for every uj £ Q, N £ N and every e G (0,1) by using Holder's inequality and the definition (52 1 of 
i?e : — > [0,oo), e G (0, 1). Therefore, the polynomial growth assumption (171 implies 



E 



< 



■ E 



L^I + \Xt\') 



+ C ■ - + R,{io) ■ N^'-'^ 



< 



iV2 



2(1 + |Xt|'* 

Me|/(^""V))|1 (i-v-«)H))-^-^ 



L • E 



. m—1 
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for every lu ^ fl, N E N and every e G (0, 1). Furthermore, we have 



E 



< 



fiXr] 



/ w 



L rn— 1 



2i • 1 + E 



\25 



for every cj e il, e N and every e G (0, 1). Using the definition (50 1 of Z: — > [0, oo) then yields 



E 



m— 1 



< 2L 

< 2L 
and finally 



\2S 



\25 



E 



< 2L 



Ar2 



1 / ^ 



1 +E 



25 



+ (C + R,{uo) + Z{uj)] .7V(^-i) 



< 2L^/CN-^ (l + E \Xt\^^ ) + (C + Re{uj) + ^(w)) • n'^'-^'^ 

< (2LVC (l + E IXtI^"^' ) + C + + • M'-^'> 



for every w G fi, G N and every e g (0, 1) due to the definition (53 1 of C. The right-hand side is finite 
according to Theorem 2.6.4 in [s]. This completes the proof of Theorem 2.1 □ 



4.4 Proof of Lemma 14.21 

Proof of Lemma Roughly speaking, tat is chosen such that the drift term alone cannot change the 
sign of the iV-th Euler approximation as long as the iV-th Euler approximation is bounded by rjv where 
N gN. Now we formalize this and observe that 

T , . T/i(0)\ ^ T . . T/i(0) 

X ■ I a; + — ■ n(x) I ^ ^ ^ — ■ ^ ■ ^^\^) ~ ^ 



N 



N 



N 



N 



>x' 



T Tm(0) 
T\x\ 



T 



a;" - ^ • ■ \^^{x) - /i(o)| 



N 
T\x\ 
N 



• \^i{x) - /i(0)| + lE\[-i,i](a;) • \^i{x) - m(0)|) 
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holds for every a; e M. Moreover, using that /i: M — > M has at most polynomial growth (see (17)) shows 



, T , . Tfi{0) 
X ■ a; H • uix) 



> x' 



N 

T\x\ 
N 



l[^i,,]{x) ■ \x\ ■ I sup^y{s)\j + lR\[_i,i](a;)-L - (2+|a;|^) j 



and hence 



> x' 



> x' 



> x' 



N 

N 

TN 
N 



• |a;| • sup |^'(s)| +lK\[_i,i](a;)-3L-|a;| 
\se[-i.i] J 



sup 



.sG[-ia] 



p^^ \^i'is)\ + 3Lj • \x\ + lR\[_i,i](a:) • 

sup |^'(s)| +3L I (\x\ + \x 
^*e[-i,il / ^ 



for every a; e M. This implies that 



N 



>x'- 



T 

N 



i sup |m'(s)|+3l) (x^ + \x\ 

Ve[-i,i] / ^ 



X^ 1 



>x' 1- 



T 
N 

T 
N 



sup |^'(s)|+3L ) (l + \x 



(5-1) 



sup \fi'{s)\+3L] {l + {rNf-'^} 
,se[-i,i] J ^ , 



holds for every x £ [— rjv,r7v]. Therefore, we finally obtain 
>x^il~^ i sup |/x'(s)| +3L ) <( 1 + max I 0, 



TV 



x^ I 1 



N 

T 

N 



,se[-i,i] 



sup |/i'(s)| + 3L 



T (sup,g[_i,i] Wis)\+3L 



N 



sel-ui] J r(sup,e[_i,i]|/i'(.)|+3i:j 



= (1 - 1) = 0. 
for every x £ [— rjv, rN]- Hence, we have shown that 



holds for all x e [0, r^r] and that 



T , , T^(0) „ 



(56) 



(57) 



(58) 



(59) 



holds for all x £ [—r^, 0). The estimates ( 58 1 and ( [59| ) give us control on the effect of the drift function in 
one direction. In addition, we need to bound the growth in the other direction. The one-sided Lipschitz 
continuity of fi (see ( [T8| ) implies 

X ■ {p{x) - /i(0)) <L-x^ 
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and hence 



X ■ {n{x) - L-x- fi{0)) < 



for every x € M. Therefore, we obtain 



for every x € [0, oo) and 



pi{x) - L • X - /i(0) < 



^(x) - L-x- /i(0) > 



(60) 



(61) 



for every x G (— oo,0). With these inequalities at hand, we now estabhsh (45) by induction on n G 
{0, 1, . . . , N} where iV G N is fixed. First of aU, we have T(j^(w) = and therefore 

\Y,^{u;)\ = \au^)\ < T 1^(0)1 + |a(0)| + |e(^)| + 1 
= e 



(Erio-fM) (T 1^(0)1 + k(0)| + |C(u;)| + 1) 

(r|M(0)| + |a(0)| + |e(c.)| + l) 



-1 



+ Y: Bgn(r,^(.))e(5:ri..."rM);3f(c.) = Dy(^).„(.) 



for every uj G flw.o, which shows (45 1 in the base case n — 0. Suppose now that (45) holds for one fixed 
n G {0, 1, 2, . . . , iV — 1}. Moreover, we fix an arbitrary uj G fl]\r,n+i C ^N,n and we now show (45) for uj 
and n + 1. For this we distinguish between the following four different cases. 
1.) First of aU suppose that (lu) > and that F„^i(cj) > holds. Then T,f (w) = T^i(a;) and 

= (^) + I • M (Y^'iu^)) + o (r„^(^)) . {W,.^^{u:) - W,.{uj)) 

T TT 

= Y^{oo) + ^ . {Y^{oo)) - L ■ y„^(c.) - ^(0)) + ^ . Y:^{u) 



and hence 



TL 



< F^Yi (a.) <Y^\u^) + — . Y^' (u) + 



N 



due to (60). Therefore, using 1 + x < for all x G M yields 

< Y^iA^) < y^{lo) .{i + ^ + a (r„^(c.)) • {w,.Ju:) - VF,«(c.))) 



= y„^(..).e""("'+^„^(a.) 

= \Y^{u:)\ ■ e"" + sgn(r„^(a;)) • /3^(..). 
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The induction hypothesis then imphes 



< Y,f_,,{u;) < <«(^),„(c.) . e"" + sgn(y„^(a;)) • /3^(c.) 
= e(^'W(^) (T 1^^(0)1 + |a(0)| + \^{u:)\ + 1) 

ri-l 

+ ^ sgn(y,^(a.)) e(^'-+^ + sgn(i;f (c)) • (3^ (u;) 



and hence 



Ziiu^)\ < e(^'W(^) (T 1^(0)1 + |a(0)| + + 1) 



(62) 



+ sgn(r,^^H)e(^"-+^""(-))/3fH 

fe=T«(w) 

= D^j^\ ^ , , (lu) — D^f} , , , (lu), 

T™{w),n+1\ ' T™_^j(Lj),ri+l ^ 

which shows that (45) holds in this case for uj and n + 1. 

2.) Suppose now that Y^{lu) < and that Y^^j^{uj) < holds. Then we also have t^{oj) T,f^i(t^) and 

T TT 

= i;f + ^ . (r„~(^)) - i • - a*(o)) + ^ • y„^(..) 



and hence 



> Y,U^) > + ^ . + ^ 



+ {<j{Y,^{oj)) - a(0)) • {w^nJco) - W,.{uj)) + a{0) ■ {w,.^^{oj) - Wt.{u;) 
due to (61). Therefore, we obtain 

+ +^(0) . {W,.ju) - w^,«h)) 

+ +^(0) . {W,.ju;) - W^,«H)) 

= - (|y„^(..)| . e"" + sgn(y„^(a.)) • /3^(a.)) . 
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Hence, the induction hypothesis yields 



= e(^'W(^) (r|MO)| + |a(o)| + ICHI + 1) 

n-l 

+ sgn(lf (c)) e(^"-+i""(-))/3f M + sgn(y„^(c.)) • ^^(w) 



5: sgn(r,^(.))e(^-..."rM)^A'(,) 
N,l I. .^ 

I 



T™(w).n+1\ / t™j(lj), 



which shows that ( 45 ) also holds in this case for ui and n+ 1. 



3.) In the third case assume that Y,^ (w) > and that F,™ ^(w) < holds. Then we obtain t„_|_]^ (w) = n+1 
Additionally, note that 

(63) 



\Y^^u;)\< sup D^fiu) 

fc,ie{0,l,...,n} 



< ''AT 



holds due to the induction hypothesis and since oj e ^N,n+i C f^Ar.n- Hence, (58) yields 



Tm(0) 
TV 



> 



(64) 



and therefore 



Y 



> 



7X0) 



and 



0>r,^i(c.)>- 



Tm(0) 



-|a(y,f(w))|- W,.^^{i^)-W,.{u,) 



> 



r|M(o)| - (L • Ir^f'l^)! + ifT(o)l) • w^t«^^(u;) - w,.{^) 



This implies 



> r„^i(w) > -T|M(0)|-(i • </(„),„(c.) + |a(0)|) •|w^,«^^(a.) - W^,«(..) 
> -T |a*(0)| - (L • tat + |a(0)|) • W^N {lo) ~ W,. (c^) 



/ Ai \ 

>-T|/i(0)|- i. — + |a(0)| • M^,«^^(u;)-PF,«(o.) 



and 



> r„^i(^) > -T 1^(0)1 - (|a(0)| + ivi) . = -T |/i(0)| - - 1 > -T |M0)| - |a(0)| - 1 (66) 

\ / TV-i A** 
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and finally 

'n+l\ 

n 

e(i:r.„+i"rM) (T 1^(0)1 + |a(0)| + 1 + leHI) + ^ sgn(y,^H)e(^r..+i"rM)^A^(^) 



|y„^iH|<TK0)| + K0)| + i + |eHI 



fc=ri+l 



which shows that (45) also holds in this case for uj and n + 1. 

4.) In the last case assume that Y^{uj) < and that Y^_^^{uj) > holds. Then we also obtain t^_^i{uj) 
n + 1. Note that 



|r^^M|< sup D^-^ic) 

fc,ie{0,l,...,n} 



holds due to the induction hypothesis and since uj £ flj\i^n+i C rjjv.n- Therefore, (59) implies 
and hence 



(67) 



(68) 



(69) 



< 



^^(0) 
N 



and 



N 



+ |a(y„^(c.))|. Wt^^(u;)-Wt^iio) 



< r|/i(0)| + {L ■ \Y^{u;)\ + 1^(0)1) • W,.Ju;) - VF,«(c.) 

< T 1^(0)1 + (l • <«(^),„(w) + |a(0)|) • \w,.Jcj) - Wt.{u;) 



(70) 



This shows 



< T 1^^(0)1 + |^L.^ + |a(0)|j •|M^,«^^(a;)-VF,«(L.) 

< T 1^(0)1 + (|a(0)| +Ni)-^=T |M0)| + ^ + 1 

and finally 

|r„^i(^)|<r|M0)| + k(0)| + i + |^(..)| 

^ e(^r.„+."rM) 1^(0)1 + 1^(0)1 + 1 + |e(c.)|) + sgn(y,^(..))e(^"-+i""(-))/3f (..) 



(71) 



fe=n+l 

(TK0)| + K0)| + i + |e(c.)|) 



+ ^ sgn(n^(..))e(^;W^""("))/3f(c.)^I?5^^(^^„^^(a.), 
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which shows that (45) holds in this case for w and n + l and which finally yields (45) for every lu e ^N,ni 
n e {0, 1, 2, . . . , TV} and every iV g N by induction. □ 

4.5 Proof of Lemma 14.31 

In order to bound the moments of the dominating process, we need to estimate the absolute moments of 
a normally distributed random variable. 

Lemma 4.7. Let il — > M fee a normally distributed T/B{R) -measurable mapping. Then we obtain 
that 

\\Yh,<p\\Y\\^, (72) 

holds for every p € [1, oo). 



Proof of Lemma \47f\ First of all, Holder's inequality implies 

\\Yh,<\\Y\\^.<p\\Y\\^,_ (73) 
for every p E [1,2), which shows (72 1 in the case p E [1,2). Denote now the mean of F: — > M by 
c := E [y] e M and the standard deviation by a := t^E (Y — cf' G [0,cx)). liY : — M is a standard 
normally distributed J^/S(]R)-measurable mapping, then 























Y 


LP 


cry + C 


< 

LP 


aY + c 


= (E 




ctY + c 




) 



< E 



E 

k=0 



\p^ 
k 



Y 



E 

> fc=0 



\p^ 
k 



Y 



^{\p^-k) 



(74) 



holds for every p £ [2,oo). Using Ely]*^ < (fc - 1)'"^^ for all fc G {2,3,...}, E|y| < ^Je{Y)^ = 1 and 
(|a| + |c|)2 < 2(a2 + c2) yields 



i^iLp < I E 



, fc=0 



\p^ 
k 



[pi 



V(bl-1)(H + N) < V^(kl + \c\) < pV{a^ + c^)^p\\Y\\ 



L2 



for every p G [2,oo), which finally shows (72). 



(75) 



□ 



Lemma 4.8. Let y ; J7 — > M be a standard normally distributed T/B{M.) -measurable mapping. Then we 
obtain that 

||e=^-l||^, <|c|e(^'+i)^' 
holds for every c G M and every p G [1, oo). 



(76) 



Proof of Lemma \4.8\ We establish (76) in the case c G (0, oo) since the case c = is trivial and since 
the case c G (— oo,0) immediately follows from the case c G (0, oo). In order to show (76 1 in the case 
c G (0, oo), note that 



E 



Ie=^-l| 



E 



E 



4y<o} 



l-e- 



<cP-E [l{Y>o} ^^"^ + • E [l{y<o} l^^ll 
= cf •E[(1{^>0} e^P^ + l{y<o}) |yn 



cY\P 



< cP {E[e^''P^] + P[y < 0]) ^ Je \Y 



|2P ^ ^2c-p- 



2^2 1 

2 



Ely 



|2p 
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and the estimate ||y||^2p < ^J2p — 1 therefore shows that 



1 \ 2p 



for all c G (0, oo) and all p G [1, oo). This completes the proof of Lemma 4. 
Lemma 4.9. Let : — > M for n G {0, 1, . . . , — 1} and N ^ N be given by (35). Then we obtain 



□ 



(77) 



for all n G {0, 1, . . . , N}, iV G N and all p G [1, oo). 

Proof of Lemma \4-S[ The triangle inequality and Lemma |4.8| impliy 



TV 




e""" - 1 


< 




LP 



N 


TL 




TL 


- 1 




— e " 


+ 

LP 


e " 


LP 


TL 
N 


e 


(W,N 




- 1 


LP 



— TJ± 

I - e N 



- y N N 

for aU n G {0, 1, . . . , - 1}, iV G N and all p G [1, oo). This completes the proof of Lemma [ij] □ 

Proof of Lemma \4.3\ Let C G (0, oo) be a real number satisfying 

5 + L + <5 + T+ 1^^(0)1 + |a(0)| <C, \^iix)\ + (x)] < C {l + \xf) (78) 

for all X G M. Such a real number C < oo indeed exists since the derivative of fi is assumed to grow at 
most polynomially according to (17). Since the exponential function is convex, we obtain that 



TL 



expK]z.a(r,^).(M/,«^-M/,«) 



(79) 



is a positive submartingale in n G {0,1,..., A^} for every z G { — 1,1} and every G N. Therefore, 
Doob's inequality (see e.g. Theorem 11.2 (ii) in (sl) shows 



E 



< 



sup e 

ne{Q,l,---:N} 



E 



p-l 



p-l 



p-l 



E 



E 



h+1 h JJE 



pzla{Y^_,)-lw^N-w^N 



(80) 



due to the moment generating function E[cxp(cy)] = exp(c^/2), c G M, of the standard normally dis- 
tributed random variable Y :— (^tJJ ~ ^tiv-i) ^^"^ hence, using |(7(a;)| < L for every x G M, 



sup e 
7ie{a,i,...,N} 



LP 



< 



p-l 



E 



; 2^* JV 



(81) 



< . . . < 



I n 



g 2^ N 



1=0 



_^]^yL-T 
p-l 
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for every p S (1, oo), z e { — 1, 1} and every N ^ N. This implies 



sup e 

ne{Q,l,...,N} 



< e 



TL 



LP 



sup e 

nG{0,l,...,Af} 



(82) 



LP 



for every p e [2,oo), 2; e { — 1, 1} and every N £N. Moreover, Lemma 4.9 shows 

5 (VTe(^'^+i)f + r) 



2 y - ~ Viv 



1 64^ 



(83) 



for all n € {0, 1, . . . , — 1}, S N and all p € [2, 00). With these estimates at hand, we now bound the 
p-th absolute moment 



E 



sup \d^:w\^ 

v,w£{0,l,...,N} 



(84) 



of the dominating process for every G N and every p £ [2, 00). By definition (37) and by the triangle 
inequality we have 



sup 

v,w£{0,l,...,N} 



and, using Holder's inequality, 



< 



LP 



sup e&^-n (C^ + iei) 

.v,we{0,l,...,N} I 



LP 



sup 

v,we{0,l,...,N} 



w—1 



^sgn(rf)e(^-V."n/?f 



LP 



sup \D': 

v,we{0,l,...,N} 



N,l 



< 



LP 



sup e 

0<v<w<N 



sup 

0<u<iii<iV 



{c' + m\L-p) 



L^p 



^sgn(rf)e&Ai"r')^ 



LP 



for all e N and all p E [2, 00). Therefore, we obtain 

< sun p,(^r=;Nr) 



sup \d^:2\ 

v,w£{0,l,...,N} 



LP 



sup ev 

0<v<w<N 



L^p 



L^P 



sup I e 

0<ti<'!i;<JV 



(E-o^af 



W — 1 



k—v 



(85) 



LP 



and 



D. 



N,l\ 



sup \-^v,w 
v,wG{0,l,....N} 



< 



sup e 
we{o,i,...,N} 



LP 



L4p 



sup e 

i;e{0,l,---,A^} 



sup e(Er=o^"r) 




sup 


we{o,i,...,N} 


L^p 


0<ti<'!i;<A' 



-(Erjo^-f) 

w—1 



(86) 



L4p 



L2p 
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for all e N and all p G [2, oo). Inequality (82 ) therefore yields 



LP 



sup \D^^'^ 

(0,1,. 



4p ^4 4p x-,4 / Q , 2p ^4 



sup 

0<t)<io<A' 



ui— 1 



sup 
e{o,i,...,N} 



5:sgn(yf)e-(^-o"n/3, 

A:— I! 

^sgn(y,^)e-(^-"r)^^ 



(87) 



L2p 



fc=0 



L2p 



for all e N and all p G [2, oo). By definition of /3,^, n e {0, 1, . . . , iV - 1}, iV e N, (see ([36])) we then 
obtain 



sup \Dv:n,\ 
i,u)e{0,l,...,Ar} 



and therefore 



LP 



, fe=0 



<e^^\c^ + ||elk^.)+2ef^^ ( Y: sgn(y,^)e-(^-"r)|^(o) 



TV' 



L2p 



+ 2e«' 



sup 
u'e{o,i,...,iV} 



^ sgn(lf ) e-(^to"r)^(o) (w'^^^^ _ w,.) 



k=0 



iN.l 



sup \D^^^ 
v,w£{0,l,...,N} 



<eP^\c^ + U\\L2p) + 2C^N- 



Lp 



+ 2Ce5' 



sup 
t«e{o,i,...,Af} 



\fe=o ^ "J 



^sgn(y,^)e-(^-"n 



fe=0 



L2p 



for all € N and all p e [2, oo). The triangle inequality and again estimate (82) hence yield 



sup l-D^ 

v,we{0,l,...,N} 



+ 2Ce5 



N,l\ 



LP 



sup 
oe{o,i,...,JV} 



W— L 

5:sgn(y-)e-(^-o^"n (e-<-l) 



fc=0 



L2p 



2Cet' 



sup 

oe{0,l,...,Ar} 



to — 1 



k=0 



L^p 



and Doob's inequality (see, e.g., Theorem 11.2 (ii) in |8 ) and Davis-Burkholder-Gundy's inequality (44) 
then show 



sup \Dv;^\ 

v,w£{0,l,...,N} 



LP 



/N-1 



, fc=0 



e-"" - 1 ) ( WtN^^ - W^ti- 



L^p 



L2p 
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for all e N and all p £ [2, oo). Holder's inequality thus gives 



sup \D^:w\ 

v,w£{0,l,...,N} 



LP 



^ fc=0 



L6p 



e-"" - 1 



P 



Wi 



Wt 



L<^p I 
N-1 



, k=0 



L4p 



and inequality (82), inequality (83) and Lemma 4.7 finally yield 

/N-l 



sup \d^:^\ 

v,w£{Q,l,...,N} 



<e'^^\c' + Uh^p) + ei^UY: 



LP 



, k=0 



6p ^4 6p ^4 1 / 



1 I X — > 8p^4 



,fc=0 



for all e N and all p E [2, oo). This shows the assertion in the case p e [2, oo). The case p E [1, 2) then 
follows from Jensen's inequality and this completes the proof of Lemma |4.3| □ 



4.6 Proof of Lemma 14.51 

Proof of Lemma \4^ First of all, we have 



sup 
ne{o,i,...,Ar-i} 



WfN^^ (lu) - WfN (oj) 



and 



and 



N-l 
n=0 



sup 

n6{0,l,...,Af-l} 



<N-P 
= 2N -1 



W.N - Wt 

Wt > N^i 



- WtN 

> N-i 
= 2N 



— Wt > N~i\ — 



sup 

2e{o,i.,...,N-i} 



Wt 



> N- 



< 2N 



2N^VT 



° 1 „zi f°° 1 x\/T _xi 

1 —^=e 2 dx < 2N I —— —e 2 dx 

K± v27r Jk± V27r iVi 

•/T VT 

xer~ dx = 2NiVT 



K± \/2n 

Vt 



s/T 



2N-WT 



2-K 



3 I \/iV 3 / -/N 



(89) 



(90) 
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for every N E N. Additionally, note that 



sup |i:'^;](a;)| > rjv > U < sup 

v,we{0,l,...,N} ' I I ne{0,l,...,Af-l} 



WjN^^ — WtN 



>N--^\ (91) 



for every N £ N. Therefore, inequality (|90| implies 



< 



sup \D^w \ > "TN 



v,w£{0,l,...,N} 



+ : 



sup 
ie{04,---,JV-i} 



3 r— \/]V 

for every e N and every p G [1, oo). Now we apply Markov's inequality to obtain 



< 



1+ sup \D^'^\\ >{l + rNY 



(92) 



E 



< 



(l + SUP^ ,^g{g 1^ JV} \Dv,W 



< sup 

\ MGN 



(l + rAr)P 

1+ sup 

u,t«e{o,i,...,A/} 



3 \/JV 

+ ^3 VTe"^ 



(93) 



LP- 



3 /— - 



and 





< sup 


H 


sup 




\ MeN 




■u,u)6{0a,...,M} 



LP/ 



< 



c 1 + sup 
\ MeN 



sup 

D,tuG{Oa,...,A/} 



> 

^-p-min(i,pi^) _^ ^|^^_^ 



LP/ 



(94) 



for every N € N and every p G [1, cx)), where we used 



^ < C - A^"'"'"(3'(7^) 

(l + rN) - 



for every A^ e N with c e (0, oo) given by c := supjygN ^Af'°"^(<i' (-s-i) )/ (i -|- ^jy)^ < oo. Moreover, (94) 
with p := 4max(4, 5 — 1) yields 

- 4max(4,(5— 1) 



< N- 



< N- 



c 1 + sup 
\ MeN 



c 1 + sup 
\ MeN 



sup \D^^'\ 

v,we{0,l,...,M} 



sup 
i>,iue{o,i,...,Af} 













ax(4,6 


-J 




i4m 


ax(4,i 







+ A^/fe"^ 



- 4 max(4,(5— 1) 



(95) 



- A^" 



T I sup a;^e 

I a;e[0,oo) 



for every N E N. The right-hand side is finite by Lemma 4.3 This proves 



d sup ( A^'' • P 
NeN 



N) 



e [0,oo). 



(96) 



Next we show that the event has probability 1. We have 



oo AI^ 



n 

M=Nm=l 



•M 



= 1 - 



oo M^ 



n n^M 







oo M^ 




)1 


= 1 -P 


U U ^^mT 


(97) 






M=Nrn=l 





24 



and hence 



n 

M=N m=l 



M 



oo 

M=N m=l 



= 1 [M 



for every N £ N due to ( 96 1 . This imphes 



oo 



> 1 - 5 



n 

M=N m=l 

for every N £ {2,3,.. .}, which shows 



ds = 1 



> 1 - 



N-l 



1 s— oo 





oo 




oo 




u 


n n^M 


= hm P 

W^oc- 


n n^M 


= 1 


Nefi 


Af=Af m=l 


Af=JV m=l 





Moreover, we have 



E 



sup 



< E 



E 



N=l 



= E 

oo 

^E 

N=l 



Ar(i+e)p 



LP 



iV(l+e)p 



(E™1i - E [In. f{Yj 



2 \ 2 ^ 



due to Lemma [4. II and therefore 



E 



sup 

NeN 



iV(l+e) 



^E 
= E 

Ar=l 



\N=1 



< {K,f [ N-^P ) ( sup E |lo„/(r^) - E I' 
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for every p g [2, oo) and every e S (0, oo). Hence, we obtain 



E 



sup 



<{2P{K,r){Y^N 

(oo 

/ oo 

< {2LK,r 



-ep 



\N=1 



sup E 



sup E 

NeN 



sup E 

NeN 



In., LP 1 



N\\P 




Y, 



N\ 
N 



s\p- 



and 



E 



sup 

NeN 



Ar(i+e) 



/ oo 

< {ALKpf J2 1 ( 1 + sup E 



\N=1 



lo„ Y^ 




< OO 



for every p € oo) and every e e (0, 1) due to Corollary 4.4 This implies 



a; e f2 



sup 

JVSN 



Ar(i+e) 



< oo 



= 1 



(100) 



for every e g (0,oo) (see also Lemma 2.1 in Kloeden fc Nc uenkirch |10| ). Putting together (98) and (100) 
shows — 1. This completes the proof of Lemma 



4.5 



□ 



4.7 Proof of Lemma 14.61 

It is somewhat inconvenient to compare the exact solution, which is a continuous time process, with the 
Euler approximations, which are time-discrete stochastic processes. Therefore, we consider the following 
interpolation process of the Euler approximation. Let Y^ : [0, T] x — > M, e N, be given by 



(101) 



for every t G n e {0, 1, . . . , TV - 1}, w e ^2 and every NeN. Note that F^^^ (cj) = Y^{lj) 

for every uj G fi, n G {0, 1, ... , N} and every NeN. Before we prove Lemma 4.6 we show that the 
restricted moments of the interpolation processes are uniformly bounded. 



Lemma 4.10. Let Y^ : [0,T] x Q. ^ R for N e N be given by ( 101 1. Then we have 



sup sup E 

NeNO<t<T 



Y 



N\P 



< OO 



(102) 



for every p £ [l,oo). 



Proof of Lem,ma \4-.10[ Inserting the definition (101 1 of the interpolation process and the polynomial 
growth of fi and a shows 



Y 



N 



< 



LP 



1+ n 



N\ 



l + LT + L\Wt~W.. 



LP 
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for every t g [t^, t^_^_i], n e {0, 1, . . . , — 1} and every iV e N. Now we apply Holder's inequality, the 
triangle inequality and Lemma |4.7| to arrive at 



N 



< 1 



LP 



< 1 



E 



N\ 



V 

^N,n I n 



L2p 
TV |2<5p' 



1 + LT + L\\Wt-Wt^ 



1+LT + L2pVf^ 



(103) 



for every t € [i^, i^+i] , {0,l,...,iV — 1} and every N G N. The right-hand side of ( 103 1 is uniformly 
bounded in n € {0, 1, . . . , N} and N G N according to Corollary 4.4 This completes the proof. □ 



Proof of Lemma Let X'^'^ : [s,T] x fl ^ R, s G [0,T], x G M, be a family of adapted stochastic 
processes with continuous sample paths given by 



xr^x+ <7ix:,ndWu 



-a.s. 



(104) 



for every t e [s, T], s e [0, T] and every a; G M. Moreover, assume that the mapping X^''{-) : R x ft ^ R 
with (x,a;) ^ X^'^^iuj) for all x e M, w e 17 is (6(]R) (g) J"t)/(S(R))-measurable for every s,t e [0,T] with 
s <t and assume that the mapping X^''{uj): M — > M with x i— >■ X^'^{uj) for all a; S M is continuous for 
every s,t G [0,T] with s < t and every oj G ft. We will show in (128) below that the difference between 

Xj"^^" and is of order 0( in a suitable weak sense for every n £ {0, 1, . . . , — 1}. Summing 

over n G {0, 1, . . . , TV — 1} for each N G N will then prove the assertion. 

First we need several preparations. According to Theorem 2.6.4 in Isl, there are real numbers Hp G 
[0, oo), p G {2, 4, 6, . . . }, such that 

E 



holds for every s,t G [0, T] with s < t and every p G {2, 4, 6, ... }, x e M. This implies 



(105) 



sup E 

t"<t<T 







p- 


= sup E 








P 










t^<t<T 













<Kp-E 



< Kp (l + E 



(i+in"r) 



(106) 



N,n n 



N\P 



for every p G {2,4,6, ...}, n G {0,1, ... , N} and every N G N. Corollary 4.4 hence implies 

Pi 



sup sup sup E 

NeN ne{0.1,....N} t^<t<T 



X. 



< OO 



(107) 



for every p G [1, oo). 

Now define u : [0, T] x M ^ M by u{t, x) = E[/(a:^'')] for every t G [0, T] and every x gR. Moreover, 
let the n-th partial derivative of u with respect to the second argument be the function u„ : [0, T] x M — M 
defined through 

Qn 



Un{t,x) 



dx^' 



u (t,x) 



for every n G {0,1, . . . ,4}, t G [0,T] and every x G M. Additionally, we use the functions u, u : [0, T] x 
R given by 

1 2 
U(t,x) = Ul{t,x) ■ /i(x) + -U2{t,x) ■ (o-(x)) 



and 



u{t, x) 



d_. 
dx 



(t,x) • ^(x) + 



1 / 92 



2 V 9x2 



U {t,x) ■ (cr(x)) 



(108) 

:M ^ 

(109) 
(110) 
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for every t E [0, T] and every a: G M. Moreover, let R E [56, oo) be a real number which satisfies 



\uit,x)\<R(^l + \xf) , 



\u2{t,x)\ <i? l + |a;| 



\ui{t,x)\<R{l + \xf^ 
\u3{t,x)\ <i?fl + |a;| 



\ui{t,x)\ <R[1 + \x\ 
\ii{t,x)\ < i? (l + |a;|^) 



\u{t,x)\ <i?(l + |a;|^ 



(111) 



dx 



u {t,x) 



<R{l + \x\ 



for every t E [0,T] and every x e M (see also Corollary 2.8.1 and Theorem 2.8.1 in |3 ). The existence of 
such a real number can be shown by exploiting (17), (18 1 and (19). In our estimates, we will need the 
real number C E [0, oo) defined by 



C := sup (N^ ■ P [(!^jv)']) + sup sup 

NeN NeNO<u<T 



N 



L + R + T + sup sup sup 

JVeNne{0,l,...,Af} t^<t<T 



X. 



<oo. (112) 



Indeed, C E [0,oo) is finite due to Lemma 4.5 Lemma 4.10 and due to (107). Moreover, since 



rN 



(113) 



holds for every t E [t„ , n G {0, 1, . . . , — 1} and every N eN, Ito's formula yields 



(.N yN \ _ f.N yN \ 



1 /■ -i+i 



u,{t^^,,Yf) {a{Y„^)f ds P- 



a.s. 



for every n E {0, 1, . . . , A'' — 1} and every N E N. Again Ito's formula yields 



+ 1 



t„+i 



+1 



t 



U2 (C+l, '^(^^f ) ^^{Y,^) dWr ds 

"3(C+i,>;^) {a{Y^))\{Y:^)drds 



1+1 



t 



1 /■ "+i 



dWr ds 



dr ds P-a.s. 



(114) 



(115) 
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for every n S {0, 1, . . . , — 1} and every N E N. 

Now we estimate all non-stochastic integrals on the right-hand side of (1151 restricted to the events 
^N.n+i for 71 e {0, 1, . . . , — 1} and N E N. For the first integral on the right hand side of (115), we 
obtain 



la 



< 



t" 



in„.„,,^^2 



LI 

dr ds 



dr ds 



and, using the polynomial growth estimate (111 I of U2, 

rti' 



la. 



+1 



< 



< 



•i+i 



"2(C+i,^;^) {^^{Y:')fdrds 
la„.„+ 



N 



1+ Y, 



N 



,25 



dr ds 



la. 



1 



Y 



N 



2+ y. 



ATI 



dr ds 



and, applying Holder's inequality and the definition (112 1 of C G [0,oo), 



+ 1 

la 

< 2L^R [ "^Y 

< 2L2i? I sup 

\ 0<u<T 



U2 



dr ds 



la. 



N 



Y 



N 



la. 



2+1: 



Af I 



ds 



L2 



1 /T 



2 V N 



for every n E {0, 1, . . . , — 1} and every e N. In addition, we have 



la. 



< 



< 



1+1 



ia„,„+.-3(C+i,i;^) (a(y„^))V(r,f) 



2L^ 



ia„.„+,u3(c+i,^;^) (i + li^n^'f) (1 



\Y, 



dr ds 



and 



la. 



< 



< 



2L3i? 



AT 



2L-'i? 



la„,„+i 1 



1; 



2+ i^. 



2+1; 



7V| 



dr ds 



dr ds 

Li 



(116) 



(117) 



(118) 



29 



and 



LI 



< 



2L^R 



< 2U'R sup 

\0<«<T 



In. 
In. 



1: 



TV 



2+ 



ATI 



N 



1 /T 

2 liV 



(119) 



for every n G {0, 1, . . . , iV — 1} and every g N. Next we use the estimates |cr(a;)| < L(l + \x\) and 
(1 + la;])" < 8(1 + x^) for every a; e M to obtain 



1 /■ "+i 



iV,,x + l 4 



< 



2L^ 



dr ds 



(120) 



LI 



and 



lnN,r. + l 



< 2L^R 



< 2L^R 



1 /■ "+i 



drds 



+1 



In. 



la„,„+i ( 1 



N 



1+ Y, 



N 



dr ds 



(121) 



Li 



N 



2+ Y 



N\ 



dr ds 



Li 



and 



la. 



< 2L^R 



1 f "+i 



+1 



dr ds 



Li 



la. 



Y 



N 



< 2L''R sup 

\ 0<u<T 



^a. 



4 DTn2/^2 i\r— 2 ^ /-r9 i\r— 2 



L2 , 



(2 + |y-r) 



dr ds 



1 /T 
2 



(122) 



< L'^RT^C^N-^ < CN- 

for every n G {0, 1, . . . , TV - 1} and every N eN. Combining \115^ , ( |118[ ), ( |119[ ) and ( |122[ ) hence yields 



E 



E 



E 



la. 



la. 



a + 1 



Ul 



U2{t::+,,Yr'')'^iyn)KY:')dWr 



ds 



1 /■ "+i 



"3(e+i,>;^) Hr,f))'dw^.ds 
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for every n G {0, 1, . . . , iV — 1} and every N £N. Due to f2jv,n+i C ^N,n we have 
and therefore 



(123) 



E[ln„,„+,w(C+i,^n+i)] -E 



W taWN , ~J,N -xrN 



N 



(^N 



< 3C"7V 



9 M-2 



E 



E 



E 



and, using that the expectation of every involved stochastic integral is equal to zero, 



+ 



f 



E 



E 



(fijV.n+l) 



+ 



2 



E 



for every n G {0, 1, . . . , — 1} and every N G^. This implies 



< 3C^Af 



'9 Ar-2 



(^^Ar,n+l) 



(f^Af,n+l) 

and, using (fijv.n+i)'^ C {SIn)" and the Ito isometry, 



1,2 



T 



< 3C^iV-2 + 



(Ojv) 



+ 



+ 



1 

IT ill ||in.,„«3(C+i,i;^) (a(y„^)) 



L2 



dr ds 



dr ds 
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for every n S {0, 1, . . . , — 1} and every N E N. Hence, ( 112 1 shows 



E[la„,„+,ii(C+i,>;+i)] 



la„„i? 1 



TV 



)i(i 



^Af I 



dr ds 



N 



+ 1 



i 1+ i: 



ds 



L2 



TV 



4L-' 1+ y, 



dr ds 



L2 



and 



E[la„,„+,u(C+i,^„+i)] -E 



RL'^CN-^ 



RLCN- 



2RL^CN-^ 



dr ds 



L2 



1+1 



AT 



2 + 



2+ 



Af I 



ds 



L2 



Af 



2+ Y 



N\ 



dr ds 



L2 



for every n G {0, 1, . . . , — 1} and every N E N. Therefore, we have 



E[lo„,„+,u(C+i,i",^i)] -E 



T 



< 3C^N~^ + RL^CN-^ 



RLCN-^ 



sup 

t" \Jt" \0<u<T 









u 


1 











sup 

tN \0<M<T 



ds 



2RL^CN-^ 



sup 

t« \Jt" \0<u<T 



In 



iv.LJW^J 



N 



dr ds 



< 3C^N-^ + RL^CN-^ / VtC ds + RLCN-^VfC + 2RL^CN-^ / VTC ds 



for every n E {0, 1, . . . , — 1} and every N E N. This finally shows that 



E[la„,„+i"(C+i,5^„+i)] -E 



< aC^iV"^ + RL^C^VtTN-^ + RLC^VtN-^ + 2RL^C^VfTN-^ 
for every n G {0, 1, . . . , — 1} and every N E N. 
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J 1 M t Y 

Next we aim at a similar estimate as ( 124) with Y^,-^ replaced by X JJ ' " for n G {0, 1, . . . , iV — 1} 
and iV e N. Ito's formula implies 



1 /■ -i+i 



M2 



ds P-a.s. 



for every n g {0, 1, . . . , — 1} and every N ^ N. This shows 



ds 



Ui 



(.N 



Xt"'^" ]a{Xs"''" ]dW. P-a.s 



rt„ .Y, 



and 



T 



■i+i 

t 



^+1 



"1 



t™ ,Y„ 



<J X. 



dWC ds 



-a.s. 



for every n S {0, 1, . . . , — 1} and every G N by Ito's formula. Hence, we obtain 



E 



< 



1 ^c.^n 



^+1 



E 



E 



E 



la 



dr ds 



E 



dx 



ds 



dW.. 



for every n € {0, 1, . . . , — 1} and every A^ e N. By (|123p we obtain 



E 



< 



E 



E 



lo„.„+i i? 1 



X 



dr ds 



1+1 



E 



E 



TV, 11 + 1 

t" 



In, 



r(X*"' " jdl^^ 



ds 



1+1 



dVK, 
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and 



E 



< RC- 



1 f T 



2 \ N 



E 



E 



9a; 



ds 



HOjv,„+i)' 



and 



E 



-E 



< -RCT^N-'^ 
- 2 



+ 



t S 



ds 

1,2 



1,2 



for every n e {0,1, . . . , N — 1} and every N gN. Therefore, we have 



and 



2 \ 2 
1,2 / 



2 , \ C- 

1,2 ; 



iV2 



+ 



and 



lnjv,n+i"(^n+l;-'^t« 



<^ + (ppivr])'M 



ln„„i? 1 



LI 



2 \ 2 

dr ] ds 

1,2 / 



E 



ln„,,„{«(C+i.5'„") + ^«(C+i,i;")} 



+ 1 



In. 



2 \ 2 

dr I ds 

1,2 



ds 
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for every n S {0, 1, . . . , — 1} and every N E N. Hence, we finally obtain 



E 



E 



< 









4 \ 2 




1 


dr ds 






J 







As 

























and 



E 



E 



dr] ds + CN-^LR 



ds 



< C^N-"^ + CN-^LR / " VTC^ds + LRC^VtN-"^ 



and hence 



E 



E 



A TV 

< C'^A^"^ + lrc^tVtn-^ + lrc^Vtn-'^ 



for every n £ {0, i, . . . , — 1} and every e N. Combining ( 124 ) and ( 127 ) yields 



E 



< lOC^A^" 



(127) 



(128) 



for every n e {0, 1, . . . , A^ — 1} and every N E N. Our interpretation of ( |128| ) is tliat the difference 
between the Euler approximation and the exact solution after a time of order 0{^) is in a weak sense 
of order 0{^). Now we split up the interval [0,T] into N G N subintervals and sum up all differences 
which arise in the subintervals. Rewriting the weak difference between the exact solution and the Euler 
approximation by a telescope sum yields 



E 



E [in 
E \ln 



E 



fixr- 



-E 



fix 



T,Y') 



■E 



-E 



E 



(129) 



and hence 



E 



E 



n = 



J2 E[lf,„./(X^"'^") 



-E 



(130) 



for every A^ e N. Since f^Ar C Q,M,n+ii we have 



(131) 
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for every n G {0, 1, . . . , iV — 1} and every N £N and therefore 



and hence 



and 



E 



■E 



N-1 



•-^ _1_1 5 ^ 



= 5^E ln„ /(X^" -/(x 

n=0 L I V 

N-1 r . , 

5]E ln„,„^, /(X^"-''")-/ X 

n — n L V V 



n = 
W-1 



n = li ^ L \ /J. 



E 



ln.-/(XT) 



-E 



E 



n = 
n = 



la.-/(XT) 



la 



N,„+iJL(aN)« 



/(x^"'^")-/[x 



r-''^ _1_1 ? J _1_ 



■E 



n = 

Af-1 



E 



■E 

n = 



la„,„+iE 



la„ • / (y;^) 

/(x^"'^")-/(x 



la. 



/(x^-^'^]-/(x;"+^ 



for every N Therefore, we obtain 



E 



< 



la.-/(XT) 

N-1 

E 



E 



n = 
N-1 



E 



N 

n + l 



la„ • / (Y^) 
( 

+ XCiV-2(||ln;..„+,-/(xf'^"') 



L2 



/ X 



n = 



L2 



la«,„+i ■ / P^: 
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and, using the Markov property and inequality (128), 



E 



N-l 



E 



< 



ilN ■ f (^W ) 



n = 

N-l 



ri = 



L2 



L2 



7V-1 JV-1 / 

n = n = \ 



J\rp 



L2 



L2 , 



for every e N. Finally, using fiAr,„+i C VlN,n for all n e {0, 1, . . . , - 1}, G N and |a::|'' <l + \x\ 
for all a; e M, we arrive at 



E 



Af-l 



E 



E 



X. 



2 AT-1 



L2 



for every e N due to ( 112 1. This proves the assertion of Lemma 4.6 
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